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Collider searches for energetic particles recoiling against missing transverse energy (^ 7 -) allow 
to place strong bounds on the interactions between dark matter (DM) and standard model (SM) 
particles. In this article we update and extend LHC constraints on effective dimension-7 opera¬ 
tors involving DM and electroweak gauge bosons. A concise comparison of the sensitivity of the 
mono-photon, mono-lU, mono-Z, mono-lU/Z, invisible Higgs-boson decays in the vector boson fu¬ 
sion (VBF) mode and the mono-jet channel is presented. Depending on the parameter choices, either 
the mono-photon or the mono-jet data provide the most stringent bounds at the moment. We fur¬ 
thermore explore the potential of improving the current 8 TeV limits at 14 TeV. Future strategies 
capable of disentangling the effects of the different effective operators involving electroweak gauge 
bosons are discussed as well. 


I. INTRODUCTION 

Because of their potential connection to DM, searches 
for represent one of the main lines of LHC research. 
These searches can be categorised based on the type of 
SM particles that recoil against DM. By now, ATLAS 
and CMS have considered a plethora of different final 
states in DM searches containing jets of hadrons, gauge 
bosons, heavy quarks and even the Higgs boson (see 
e g- [D for a recent review of the experimental status). 

In most cases these studies are performed in the con¬ 
text of an effective field theory (EFT) which correctly 
captures the physics of heavy particles mediating the in¬ 
teractions between DM and SM fields, if the mediators 
are heavy enough to be integrated out. Below we will 
consider the effective Lagrangian 

E W 

k=B,W,B,W 

which contains the following four SU{2)lxU(1)y gauge- 
invariant dimension-7 operators 

Ob = XX , Ow^ = XX , 

( 2 ) 

Os = XXB^uB^\ = xxWl.W^'^^''. 

Here A represents the scale of new physics at which the 
higher-dimensional operators 0 are generated, i.e. the 
scale where the messenger particles are removed as ac¬ 
tive degrees of freedom. The DM particle y can be 
both a Dirac or a Majorana fermion and — 

d.B^ iK- = - d^Wl + g 2 €^^^W/,Wl!) IS the 

U{1)y (S'?7(2)i) field strength tensor, while = 

^/‘^'efj.uXpB^P = l/2e^^ApW^*’^'’) denotes its dual 

and g 2 is the weak coupling constant. 
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FIG. 1: Representative examples of graphs that generate a 
4 “ T 5 Z (— y t -\- WjZ (—'f 7 ) or a + 2f sig¬ 

nal. The operator insertions are indicated by yellow squares 
while SM vertices are represented by black dots. Propagators 
labeled by V include all possible photon, Z-boson or VE-boson 
exchanges. See text for further details. 


The operators introduced in 0 appear in models of 
Rayleigh DM (see for instance [SB]). They are special 
in the sense that, up to dimension 7, they are the only ef¬ 
fective interactions which lead to velocity-suppressed an¬ 
nihilation rates of DM to photon pairs While the 

sensitivity of future direct detection experiments may al¬ 
low to set novel bounds on the Wilson coefficients Cb(A) 
and Cw (A) for heavy DM particles with ^ 1 TeV 
once loop effects are taken into account [8], in the case 
of light DM the leading (and for C^{A) and C^{A) the 
only) restrictions arise and will continue to arise from 
collider searches involving large amounts of f^ct, 

the DM-SM interactions 0 have been constrained us¬ 
ing 7 TeV and 8 TeV LHC data on invisible decays of the 
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Higgs boson in the VBF mode [9] as well as the ^ 

[TOl [11] , the mono-photon [12] and the + W [IS] chan¬ 
nel. 

The main goal of this article is twofold. First, to up¬ 
date the existing constraints by taking into account the 
latest mono-photon [T4l[T5] . i]rj.-\-W/Z hadrons) [T6] 
and VBF h invisible m searches. Second, to ex¬ 
tend the studies leHiS] by considering in addition the 
leptons) channel [T8j[T9] as well as the newest 
mono-jet data [20]. An assortment of Feynman diagrams 
that lead to the signatures investigated in the follow¬ 
ing are displayed in Fig. [^ 

Our analysis shows that depending on the choice of pa¬ 
rameters, either the mono-photon or the mono-jet data 
give rise to the strongest restrictions at present. By com¬ 
bining the information on all available channels we are 
thus able to derive bounds on the coefficients Ck{A.)/A^ 
in 0 that improve on the existing limits. Building 
upon [21], we furthermore demonstrate that measure¬ 
ments of the jet-jet azimuthal angle difference in $j^-\-2j 
events may be used to disentangle whether the DM bilin¬ 
ear XX couples more strongly to the combination 

or the product of field 

strength tensors. Similar ideas have also been brought 
forward in [9]. 

The outline of this article is as follows. In Sec. nn we 
review the existing LHC searches for signatures that 
we will use to constrain the effective interactions 0. 
In Sec. |HI| we derive the restrictions on the parameter 
space by combining all individual search modes, com¬ 
menting also on how future measurements may improve 
these limits. This section contains in addition a discus¬ 
sion of the azimuthal angle correlations between the two 
jets in the $rj. + 2j channel. Our conclusions are pre¬ 
sented in Sec. IlYl 


ferent signal regions with a varying cut on the trans¬ 
verse momentum of the photon (pt, 7 )- Note that due 
to the higher-dimensional nature of the operators 
the $rj. + 7 signal has a rather hard pt ,7 spectrum. As a 
result, we find that the most severe cut ofpT ,7 > 700 GeV 
gives the strongest bounds on the parameter space in our 
case. The corresponding 95% confidence level (CL) limit 
on the fiducial cross section reads 

0’fid(pp $T+l) < 0-22 fb . (4) 


B. Mono-Z signal 


In the case of the $rj. -\- Z search channel, we 

use the ATLAS results [H], that utilise 20.3 fb“^ of 8 TeV 
data. The selection criteria relevant to our analysis are 


Pt/ > 20 GeV , \r]i\< 2.5 , mu G [76,106] GeV, 

1^ I ^ o c: \PTM - $t\ ^ n c: 

2.5 , < 0.5 . 

Pt,u 


Here mu^ Tju and Pt,u denote the invariant mass, the 
pseudorapidity and the transverse momentum of the di¬ 
lepton system, respectively. The ATLAS analysis defines 
four signal regions with different lower thresholds. As 
it turns out, in the considered case the requirement > 
350 GeV gives rise to the best bounds. Including Z-boson 
decays to both electrons and muons {i = e,/i), the AT¬ 
LAS experiment obtains for this cut the following 
95% CL bound 


cTfid {pp —^ T ^ ^ ^ )) 0.27 fb . (6) 


C. Mono-LF signal 


II. SEARCH CHANNELS 

In this section we list the various cuts and the values 
of the fiducial cross section (crfid) of each individual $rj. 
channel. This information will be used in the next section 
to set limits on the coefficients Ck{A)/h? appearing in 
the effective Lagrangian 0 . 

A. Mono-photon signal 

We begin with the mono-photon signal, which has re¬ 
cently been searched for by both CMS [14] and AT¬ 
LAS [15]. Since the former search leads to the stronger 
restrictions, we employ the CMS results, which are based 
on 19.6 fb“^ of 8 TeV data. The relevant cuts are 

$T > 140 GeV, \r]^\< 1.4442 , (3) 

where denotes the pseudorapidity of the photon. The 
CMS collaboration performs the measurement in six dif- 


Both ATLAS [18] and CMS [19] have searched for a 
mono-VF signal in the leptonic decay mode. We find that 
the ATLAS search for the final state, which uses 
20.3 fb“^ of 8 TeV data, gives the strongest constraints, 
and thus we consider only this channel. The most im¬ 
portant experimental cuts are 

PT,M>45GeV, \r]^\ G [0,1] U [1.3, 2], 

/- 

rriT = ^2pT,n$T (l - , 

where mr is the transverse mass which depends on the 
angle between the pt^^ and the vectors. AT¬ 

LAS sets bounds on cTfid for three different ttit cuts, and 
like in the case of the mono-photon signal, we observe 
that the strongest restriction of mr > 843 GeV provides 
the best limits on the interactions 0 . At 95% CL the 
bound on the corresponding fiducial signal cross section 
is given by 

CTfid {pp $T + W pvp) < 0.54 fb. (8) 
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D. Mono-W/Z signal 

The ATLAS search [16] looks for a $rj. jz signal, 
where the W or Z boson decays hadronically. This anal¬ 
ysis is based on 20.3 fb~^ of 8 TeV data, jet candidates 
are reconstructed using the Cambridge/Aachen (C/A) 
algorithm [22] with a radius parameter R= 1.2 and sub¬ 
jected to a mass-drop filtering procedure [23]. Events are 
required to have at least one C/A jet with 

PTj > 250GeV, \r]j\ < 1 . 2 , 

(9) 

rtij e [50,120] GeV , ^ > 0.4. 

Here rrij refers to the mass of the large-radius jet, while 

\/y = min(pr,ii,PT,iJ is a 

measure of the momentum balance of the two leading 
subjets ji and j 2 contained in the C/A jet. The 95% CL 
limits on the fiducial cross section depend also on the 
imposed threshold, and it turns out that the stronger 
of the two cuts, i.e. > 500 GeV, provides the most 
stringent constraints. In this case, the relevant limit on 
the fiducial cross section is 

{pp $T + W/Z hadrons)) < 2.2 fb . (10) 

E. Mono-jet signal 

One can also use mono-jet events to constrain the op¬ 
erators in since the corresponding searches allow for 
the presence of a secondary jet. Here we will employ 
the newest GMS results m, which make use of 19.7 fb ^ 
of 8 TeV data. Like GMS, we reconstruct jets using an 
anti-Zct algorithm [24] with radius parameter R = 0.5. 
The relevant selection cuts are 

PTJ, >110GeV, <2.4, 

PT,h > 30 GeV, |?7 j2 | < 4.5, (n) 

V 2.5, 

where ^(t)j,j2 is the azimuthal separation of the two lead¬ 
ing jets. Another important selection criterion is the 
imposed jet-veto [26], which rejects events if they con¬ 
tain a tertiary jet with Ptj^ > 30 GeV and 177 ^ 3 ! < 4.5. 
The GMS measurement is performed for seven differ¬ 
ent $T regions, and we find that for the considered in¬ 
teractions the highest sensitivity is obtained for $t > 
500 GeV. The corresponding 95% GL limit on the fidu¬ 
cial cross section reads 

r^fid(pp ^ $T + 2jf) < 6.1 fb . (12) 

F. VBF invisible Higgs-boson decays 

Last but not least, we consider the results of the CMS 
search for invisible decays of the Higgs boson in the VBF 


channel which uses a 8 TeV data sample, corre¬ 

sponding to an integrated luminosity of 19.5 fb“^. Jets 
are reconstructed employing an anti-/ct clustering algo¬ 
rithm with R = 0.5, and subject to the following require¬ 
ments 

PT,h > PT,h > 50 GeV, \r]j^ I, his I < 4-7, 

Vji ■ Vj2 ^ 0, Arij-^j2 > 4.2, (13) 

rrij^j^ > 1100GeV, ^4^jij2 < 1-0 • 

The missing-energy cut is $t >130 GeV and a central 
jet-veto is imposed to any event that has a third jet with 
PtJs > 30 GeV and a pseudorapidity between those of 
the two tagging jets. For these cuts, CMS obtains the 
following 95% CL bound on the fiducial cross section 

cTfidiPP ^ 2j) < 6.5fb. (14) 


III. NUMERICAL RESULTS 

In order to determine the cross section for the $t sig¬ 
nals associated to the effective operators (^, we have 
implemented each of them in FeynRules [27 , generating 
a UFO output [28] . The actual event generation has been 
performed at leading order with MadGraph 5 [29] utilis¬ 
ing CTEQ6L1 parton distributions [30]. Parton-shower 
effects and hadronisation corrections have been included 
by means of PYTHIA 8 [31] and jets constructed us¬ 
ing Fast Jet 3 [32]. We employ Delphes 3 [25] as a fast de¬ 
tector simulation to estimate the reconstruction efficien¬ 
cies for the different '$t signals. The efficiencies that we 
find amount to around 70% for the mono-photon signal, 
60% in both the mono-Z and mono-IT case and 65% for 
the mono-kF/Z signature. These findings agree with [15] 
for the $iT + >5 [H] [O] 

and m for the $T -PW/Z signal. For the mono-jet sig¬ 
nal and the search for invisible decays of the Higgs boson 
in the VBF channel, we find reconstruction efficiencies in 
the ballpark of 95%. 

Our Monte Garlo (MG) implementation has been val¬ 
idated by reproducing the numerical results of mm 
within theoretical uncertainties. These errors have been 
assessed by studying the scale ambiguities of our re¬ 
sults. We have used the default dynamical scale choice 
of MadGraph 5, varying the scale factor in the range 
[ 1 / 2 , 2 ]. We find that the predictions for the mono¬ 
photon, $T + ^ (^ i^i~) and $T + ^ (^ cross 

sections calculated in this way vary by around ±15%, 
while in the case of the $t + W/Z (^ hadrons), the 
mono-jet and the VBF h invisible signal, relative dif¬ 
ferences of about ±20% are obtained. Note that these 
errors are smaller than those found in [26] [33H35] , since 
all the tree-level $t cross sections considered in our work 
do not explicitly depend on as . The quoted uncertainties 
thus reflect only the ambiguities related to the change of 
factorisation scale, but not renormalisation scale. 
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FIG. 2: Assortment of LHC bounds on the new-physics 
scale A, assuming Cs(A) = 1, Cw{^) = 0 (upper panel) 
and Cs(A) = 0, Cu^(A) = 1 (lower panel). In both cases the 
DM particles are taken to be Dirac and G^(A) = C^(A) = 
0. The coloured curves correspond to the limits arising from 
the latest mono-photon (red), + Z l~) (orange), 

$rj. + W (^ (yellow), $rp -^WIZ (^ hadrons) (purple), 
mono-jet (blue) and VBF h invisible (grey) searches. The 
width of the bands reflect the associated scale uncertainties. 


A. Dependence on a single Wilson coefficient 

In Fig. we present the limits on the new-physics 
scale A for C^{A) = C^{A) = 0 and the two choices 
Cb{A) = 1 , Cw{A) = 0 (upper panel) and Cb{A) = 0 , 
Cw{A) = 1 (lower panel) for the Wilson coefficients eval¬ 


uated at A. The shown predictions correspond to Dirac 
DM and the widths of the coloured bands illustrate the 
impact of scale variations. For Cb{A) = 1 , Cw{A) = 0 , 
one observes that the mono-photon search la provides 
the strongest constraints in most of the parameter space. 
Numerically, we find that the scale A has to satisfy A > 
510 GeV for ^ 100 GeV in order to meet the 95% CL 
requirement (3. In the case Cb{A) = 0 , Cw{A) = 1, 
on the other hand, the latest mono-jet data [ 20 ] impose 
the leading restrictions. At 95% CL, the inequality 
translates into a lower limit of A ^ 600 GeV for DM 
masses below 100 GeV. The shown limits also hold in 
the case that C^{A) = 1 , C^{A) = 0 or C^{A) = 0 , 
C^{A) = 1 and Cb{A) = Cw{A) = 0 , while for 
Majorana DM the constraints on A would be stronger 
by around 12 %. Note finally that ^ 

searches do not provide any constraint on scenarios with 
Cw{A) = C^{A) = 0 . 

To better understand the restrictions imposed by the 
various search channels, we consider the Feynman rules 
associated to the effective operators Ob and Ow enter¬ 
ing 0 - In momentum space, the resulting interactions 
between pairs of DM particles and SM gauge bosons take 
the form 

4? 

9ViV2 {pTp2^ - Pi • P 2 ), (15) 

where pi {pi) denotes the momentum (Lorentz index) of 
the vector field Vi and for simplicity the spinors associ¬ 
ated with the DM fields have been dropped. In terms 
of the sine (s^^;) and cosine (c^j) of the weak mixing an¬ 
gle and the Wilson coefficients Cb{A) and (7w(A), the 
couplings gViVj read 

gAA = c^Cb(A) + sl^CwiA) , 

Qaz = -SwCw {Cb{A) - Cw{A)) , 

(16) 

gzz = 5^Cb(A) + c^Cw{A ), 
gww = Cw{A). 

These results do not coincide with the expressions re¬ 
ported in unmans]- From we see that in the cou¬ 
pling qaa of DM to two photons, the Wilson coefficients 
Cb{A) enters compared to Cw{A) with a relative factor 
of (?^l 3.3. On the other hand, in the case of the 
coupling between DM and Z-boson pairs gzz^ the depen¬ 
dence on Syj and Cy^ is reversed compared to These 
properties explain why the limit on the new-physics scale 
A from mono-photon ^ Z ^WjZ 

hadrons), mono-jet and VBF h ^ invisible) searches is 
stronger (weaker) in the upper panel than in the lower 
panel of Fig. 

A second important feature worth noting is that chan¬ 
nels with leptons in the final state typically lead to weaker 
restrictions on the parameter space than modes involving 
hadrons. This is a simple consequence of the fact that the 
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Cb(A) Cb(A) Cb{A) 

FIG. 3: Limits on A in the Cb{A)-Cw{A) plane. The different panels correspond to the mono-photon (upper left), (!rp-\-Z 

(upper middle), + W (^ (upper right), + WjZ (^ hadrons) (lower left), mono-jet (lower middle) and VBF 
h invisible (lower right) search. All results employ — 100GeV and C^{A) = C^(A) = 0. The contour labels indicate 
the value of the new-physics scale in units of GeV. 


electroweak SM gauge bosons dominantly decay hadron- 
ically. Numerically, one has Br (Z ^ ^+^“) 7% and 

Br (LF ^ iiv^) 11%, while Br (Z ^ hadrons) 70% 
and Br {W hadrons) 68% [36]. The strong suppres¬ 
sion of the leptonic decay widths overcompensates the 
higher detection efficiencies of final states involving lep¬ 
tons, and as a result the LHC searches for + hadrons 
are superior to those looking for + leptons signals. 

Our third observation is that the latest mono-jet data 
are evidently more constraining than the recent VBF h 
invisible search. While these analyses explore the same 
final state, i.e. they probe quite different parts of 

the phase space. In fact, the selection criterion that has 
the biggest impact in our study is the rather loose missing 
transverse energy cut of > 130 GeV imposed in the 
VBF h invisible search. This selection is tailored for 
a Higgs boson of 125 GeV, but fares less well if one tries 
to probe higher-dimensional operators of the form 
Since the operators Ob {O^) and Ow {O^) produce a 
rather hard spectrum, more severe requirements 
allow for a cleaner separation between signal and SM 
background. 


B. Dependence on two Wilson coefficients 

Until now we have have studied the constraints on the 
new-physics scale A as a function of the DM mass m^, 
keeping the values of the high-scale Wilson coefficients 
fixed. In the panels of Fig. we instead show contours 
of constant A in the (7 b(A)-(7w(A) plane. In all plots 
we employ = 100GeV and set C^{A) = C^{A) = 0. 
The first noticeable feature of the shown predictions is 
that only the mono-photon signal depends more strongly 
on Cb{A) than (7w(A), while for all the other chan¬ 
nels the situation is reversed. Second, with the excep¬ 
tion of the mono-photon case, one observes that the ma¬ 
jor axes of the elliptic contours in all panels are almost 
aligned with the Cw{A) axes. This means that interfer¬ 
ence effects between contributions arising from Ob and 
Ow are small in all of these cases. The third impor¬ 
tant property following from the colour shading of the 
depicted results is that currently either the newest mono¬ 
photon or the mono-jet data provide the leading bounds 
in the entire Cb{A)-Cw{A) plane. This feature is fur¬ 
ther illustrated by the upper panel in Fig.[^ In this plot 
the overlaid numbers indicate the search strategy that 
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FIG. 4: Combination of the bounds on the new-physics scale 
in the Cs(A)-Cu^(A) plane, employing = 100GeV and 
C^{A) = C^{A) = 0. In the upper panel the search strategy 
that provides the leading constraint is indicated by the super¬ 
imposed numbers, with 1 (5) representing the latest mono¬ 
photon (mono-jet) search, while the lower panel shows the 
resulting contours of constant A in units of GeV. 
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contributes the best sensitivity on A at each point, with 
1 and 5 corresponding to the mono-photon and mono¬ 
jet channel, respectively. One sees that if the ratio of 
Wilson coefficients satisfies |(7 b(A)/(7vp(A)| ^1.5 then 
the limit @ gives rise to the strongest constraint, while 
in the remaining Cb(A)-Cvp(A) plane the bound (14) is 
most restrictive. The A contours obtained by combining 
all available channels are depicted in the lower panel 
of Fig. [§ 


Finally, we note that the values of A that can be ex¬ 
cluded with the current data are low compared to typical 
LHC energies. In order to go beyond the EFT descrip¬ 
tion, one has to specify a ultraviolet (UV) completion, 
where the operators in (|^ arise from a renormalisable 
theory after integrating out the heavy degrees of free¬ 
dom mediating the interactions. UV-complete models 
that generate the operators Ob and Ow through loops 
of states charged under U{1)y and/or SU{2)l have been 
proposed in [3] and their LHC signatures have been stud¬ 
ied in 0. If these new charged particles are light, the 
high-pT gauge bosons that participate in the Pro¬ 
cesses considered here are able to resolve the substructure 
of the loops. This generically suppresses the cross sec¬ 
tions compared to the EFT predictions [33] , and thus will 
weaken the bounds on the interaction strengths of DM 
and the electroweak gauge bosons to some extent. Fur¬ 
thermore, the light charged mediators may be produced 
on-shell in pp collisions, rendering direct LHC searches 
potentially more restrictive than searches. Making 
the above statements precise would require a study of a 
concrete UV completion. 


C. Future sensitivity 


It is also interesting to explore how the reach on the 
new-physics scale A might improve at the 14TeV LHC. 
In what follows, we will only consider the mono-jet signal, 
applying the event selection criteria that have been used 
in the sensitivity study by ATLAS m- These read 


Ptji > 300 GeV , \r]j^ | < 2.0 , 

PT,j2>50GeV, \T]j^\<3.Q, 

> 0.5, 


(17) 


and jets are reconstructed using an anii-kt algorithm with 
R = OA. Events with a third jet of Prjs > 50 GeV and 
\r]j^\ < 3.6 are vetoed and the missing transverse en¬ 
ergy cut that we employ is > 300 GeV. Note that 
compared to (11) the PtJi^ PT,j 2 thresholds are 


increased both to avoid pile-up and to enhance the signal- 
over-background ratio. In order to determine the limits 
on the scale A, we take crfid(pp ^ ^ + j) = 

5.5 fb [37], assuming a total systematic uncertainty on 
the SM background of 5%. Eor the choice Cb{A.) = 0, 
Cw{^) = 1 and Cp{A) = C^{A) = 0, we find that with 
25 fb“^ of data, corresponding to the first year of running 
after the LHC upgrade to 14 TeV, one may be able to set 
a 95% CL bound of A > 1.3 TeV for < 100 GeV. 
Compared to the present limit, this corresponds to an 
improvement of the bound on A by more than a factor 
of 2. With 300 fb“^ and 3000 fb“^ of accumulated data, 
we obtain instead A > 1.4 TeV. These numbers make 
clear that at 14 TeV the sensitivity of + j searches 
will rather soon be limited by systematic uncertainties 
associated to the irreducible SM background. To what 
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extent this limitation can be evaded by an optimisation 
of the mono-jet searches and/or an improved understand¬ 
ing of the pp Z vv) + j channel, would require a 
dedicated study. Such an analysis is beyond the scope of 
this work. 


D. Analysis of jet-jet angular correlations 


So far we have analysed only observables that are in¬ 
sensitive to whether the $rj. signal is generated by an in¬ 
sertion of the effective operator Ob (Ow) Op {O^). 
This ambiguity can however be resolved by measuring 
the azimuthal angle difference of forward jets pro¬ 

duced in $rj. 2j events [9l|2T]. Besides the cuts (17), 
we impose the following VBF-like selection requirements 
in our analysis 


Vji ■ Vj 2 < 0 > ^Vjij 2 > 2, > 1100 GeV. (18) 


Here the cut on the pseudorapidity separation helps to 
sculpt the angular correlations between the tagging jets, 
while the di-jet invariant mass threshold improves the 
signal-over-background ratio. 

In order to understand why the operators Ob {Ow) 
and Op {O^) lead to different jet-jet angular correla¬ 
tions, one has to consider their Feynman rules. In the 
case of the operators containing regular field strength 
tensors this has already been done in (15), while for their 
dual counterparts we obtain 


^ 9ViV 2 {piuP2X - Pi\P2v) , (19) 

with gViVj given in ( [l^ . The selection cuts ( [l8| ) em¬ 
phasise the parts of the phase space where the exter¬ 
nal partons experience only a small energy loss and the 
momentum components of the tagging jets in the beam 
direction are much greater than those in the transverse 
plane. In this limit the structure of the pp 
matrix elements is straightforward to work out [38] . 
In the case of the effective operator Owi one gets for 
instance A4w ^ -92-91 MI-P 2 /X 2 ) 

PtJi 'PT,j 2 ^ while for O^ one arrives instead nt A4^ ^ 

<^t^i(i 2 i'xJi"J 2 ^PiP 2 ~ PT,ji X pT,i 2 - Here J* and Pi de- 
note the currents and momenta of the electroweak gauge 
bosons that partake in the scattering. These simple ar¬ 
guments imply that the /S.(j)j^j^ spectrum corresponding 
to Ow should be enhanced for collinear tagging jets, 
A0jij2 = 0, while for = 7r/2 it should show an 

approximate zero. In the case of O^, on the other hand, 
the ^2 distribution should have a dip if the two jets 
are collinear, = 0, or back-to-back, /^(t)j^j 2 = tt. 

Note that the above arguments do not depend on the chi¬ 
rality of the DM current. This means that Ob^ Owi Op^ 
O^ and the operators obtained from ([^ by replacing 
XX with X 75 X lead to very similar jet-jet angular corre¬ 
lations, as we have explicitly verified. 

In Fig. we plot the A(t)j^j^ spectra for the choices 
Cb{A) = 0, Cw{^) = I (red curve) and Cp{A) = 0, 



FIG. 5: Azimuthal angle distributions at the 14TeV LHC. 
The signal curves correspond to Cb{A) — 0, Cw{A) — I (red) 
and Cp{A) — 0, C^(A) = 1 (blue), and both use A = 1 TeV 
and = 100 GeV. For comparison the prediction of the 
dominant SM background process pp ^ Z (^ z/z/) + 2j (black) 
employing the same event selection criteria is shown as well. 


C^{A) = I (blue curve). All shown predictions are ob¬ 
tained for the 14 TeV LHC and employ A = I TeV and 
= 100 GeV. The fiducial signal cross sections amount 
to I.Ofb, independently of whether the insertion of Ow 
or O^ is considered. The expected sine-like (cosine-like) 
behaviour of the modulation in the azimuthal angle dis¬ 
tribution corresponding to Ow {Ow) is clearly visible 
in the figure. These shapes should be contrasted with 
the spectrum of the dominant SM background process 
pp ^ Z {-^ Pu)-\-2j (black curve), which is rather fiat for 
values A(j)j^j^ ^2.6 and then rapidly drops to zero. The 
corresponding fiducial cross section is 0.35 fb, implying a 
signal-over-background ratio of S/^/B 8.4, 29 and 93 

for 25fb“^, 300fb“^ and 3000 fb“^ of data, respectively. 

The given S/\/B values imply that running the LHC 
for a couple of years at 14 TeV should provide a sufficient 
number of events to analyse the jet-jet angular correla¬ 
tions. To quantify this statement, we use a toy MC and 
generate event samples for both signals and background 
corresponding to 300 fb~^ and 3000 fb~^ of luminosity. 
The resulting differential cross sections are then fitted 
to [39] 


I da 


2 

an cos . 

n=0 


( 20 ) 


The coefficient ao is fixed by the normalisation of the 
A(j)jij^ spectrum, and the ratio ri = aija^ turns out to 
be rather insensitive to which type of higher-dimensional 
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interactions is considered. In contrast, the combination 
^2 = ^^ 2 / 0^0 is a measure of the CP nature of the inter¬ 
actions that lead to the 2j final state (see e.g. [38l l39] b 
This ratio is expected to be positive (negative) for an 
insertion of Ob {O^) and Ow (O^). We stress that 
by considering normalised distributions, theoret¬ 

ical uncertainties are reduced and that the predictions 
become fairly independent of EFT assumptions m- 
In Fig. we present the results of our toy MC. 
The upper panel (lower panel) corresponds to 300 fb~^ 
(3000 fb“^) of LHC data collected at 14TeV. The 
expected azimuthal angle distributions for the signal 
plus background predictions are coloured blue (red) for 
Ow (O^). For comparison, the SM-only result (grey) 
divided by a factor of 3 is also shown. The solid curves 
illustrate the best fits to ( [^ , restricting the rapidity sep¬ 
aration A(j)j-^j^ to the range [0, 2.5]. For 300 fb~^ of data, 
we obtain for r 2 the central values and uncertainties 

(^ 2 )w+SM ~ ^ ’ 

(^2)^+SM ~ ^ 

(^2)sm =-0-12 ±0.22. 

In the case of 3000 fb~^ of luminosity, we find instead 


(^2)w+SM ~ ^ ’ 

(^ 2 )^+SM ~ “0.40 ± 0.04, (22) 

(^ 2 )sm =-0-13 ±0.07. 


We observe that for Ow {O^) the combination r 2 is 
indeed positive (negative). Defining a significance as 
Sk = ((^2)/c+SM - (^2)sM)/(Ar2)/c+sM, wc get from (|^ 
the values sw = 2.7 and = —2.4, while (22) leads to 
sw = 10.3 and = — 6 . 8 . Our toy MC study corre¬ 
sponding to 300fb“^ (3000 fb“^) of data hence suggest 
that a distinction between the azimuthal angle distribu¬ 
tions of Ow and O^ at the 5 cr (17cr) level should be 
possible at the 14TeV LHC. We emphasise that our toy 
study assumes a perfect detector and that we have not 
optimised the cuts (18) to achieve the best significance. 
Once the data is on tape, it will become an experimen¬ 
tal issue of how stringent the VBF-like selections can be 
made to extract the most information on the jet-jet an¬ 
gular correlations for a given limited sample size. 


IV. CONCLUSIONS 

In this article we have studied LHC constraints on ef¬ 
fective dimension-7 operators that couple DM to the SM 
electroweak gauge bosons and emphasised the comple¬ 
mentarity of different searches for constraining the 
associated Wilson coefficients. Focusing on the interac¬ 
tions that induce only velocity-suppressed annihilation 
rates, we have combined the information on all individual 
search modes that are available after LHC run- 1 . In this 



FIG. 6: Normalised distributions for 300 fb“^ (up¬ 

per panel) and 3000 fb“^ (lower panel) of 14TeV LHC data. 
The red (blue) histogram shows the signal plus background 
prediction for Ow (Ow)- The grey bar chart represents the 
expected SM background, which for better visibility, has been 
rescaled by a factor of 1/3. The solid curves indicate the best 
fits of the form ao -h ai cos A(j)j^j^ -h a 2 cos See text 

for additional explanations. 


way we are able to derive bounds on the new-physics scale 
A that exceed all previous limits. Our studies show that 
at present, depending on the choice of parameters, ei¬ 
ther mono-photon or mono-jet searches provide the most 
severe constraints on the considered dimension-7 interac¬ 
tions. For DM masses ^ 100 GeV and Wilson co¬ 
efficients |(7/c(A)| 1 , the existing 8 TeV LHC searches 
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allow to exclude values of A below about 600 GeV at 
95% CL. The improved reach of analyses in 2015 
and beyond is also studied, finding that with 25fb~^ of 
14TeV data, LHC mono-jet searches should be able to 
improve the latter bound to approximately 1.3 TeV. Be¬ 
yond this point further progress will be hindered by the 
imperfect understanding of irreducible SM backgrounds 
such 8iS pp ^ Z vv) + j. Finding ways to overcome 
these limitations will be crucial to exploit the full physics 
potential of searches to be carried out at later stages 
of the LHC. 

We have furthermore emphasised that given the large 
statistics expected at the phase-1 and phase-2 upgrades 
of the 14 TeV LHC, searches should be able to not 
only determine integrated, but also differential cross sec¬ 
tions. From the theoretical point of view, such nor¬ 
malised distributions have the clear advantage, that com¬ 
pared to the total cross sections theoretical uncertainties 
are reduced and that the obtained predictions depend 
only weakly on the assumptions underlying the EFT de¬ 
scription. As an example we have explored the prospects 
to measure jet-jet angular correlations in $rj.-\-2j events. 
Taking into account the pseudorapidity correlations of 
the two tagging jets, the resulting distributions in the 
azimuthal angle separation exhibit the relative 

strength of CP-even and CP-odd interactions of DM with 
gauge boson pairs. Our toy MC studies indicate that 
already with 300 fb“^ of data a distinction between the 


new-physics and the SM-only hypotheses can be achieved 
at a statistically significant level, and that the sensitivity 
of the discussed searches is greatly improved by going to 
3000 fb~^ of luminosity. A more precise determination of 
the analysing power, including systematic uncertainties 
would require a full detector simulation, which is beyond 
the scope of the present article. We however believe that 
it is imperative that the ATLAS and CMS collaborations 
direct some activity towards the study of differential dis¬ 
tributions of final states like + 2j. 
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